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SUMMARY 

This  report  describes  the  results  of  a  study  that  was  conducted  to  determine  the  effects  of 
key  experimental  parameters  (such  as  the  amount  of  energy  deposited,  the  panel  thickness,  the 
depth  of  the  discontinuity,  the  normal  thermal  conductance,  and  the  lateral  thermal  conductance) 
on  the  peak  thermal  contrast.  For  that  purpose,  a  panels  with  various  flat  bottom  holes  was 
fabricated.  Different  amounts  of  heat  were  deposited  on  the  surface  of  the  panel  and 
thermographic  data  were  gathered.  A  simple  model  was  developed  which  takes  into  account 
lateral  heat  conduction,  anisotropic  thermal  conductivity,  thickness,  flaw  size,  density,  and  pulse 
duration.  The  final  model  provides  a  simple,  yet  accurate,  description  of  the  time  dependence  of 
the  thermal  contrast. 


INTRODUCTION 

As  active  thermography  continues  to  gain  acceptance  as  an  NDE  tool  for  both  in-service 
and  manufacturing  applications,  the  need  for  robust  physical  models  of  the  process  increases. 
Although  the  sophistication  of  both  analytical  and  numerical  modeling  efforts  has  increased 
considerably1,2’3,  development  of  a  3-dimensional  model  for  anisotropic  media  is  still  a 
formidable  task. 

In  a  previous  paper4,  a  simple  model  was  developed  (referred  to  as  Model  I)  to  describe 
the  surface  temperature  evolution  of  a  1/8  in.  thick  aluminum  panel  with  fixed  1  in.  diameter  flat 
bottom  holes  of  different  depth  in  response  to  a  short  pulse  of  radiant  energy.  That  model 
correctly  described  the  temporal  behavior  of  all  1  in.  flat  bottom  holes.  To  further  validate  Model 
I,  a  new  panel  was  fabricated  (figure  1)  that  contained  flat  bottom  holes  of  different  diameters 
and  different  depths.  After  careful  analysis  of  the  data,  it  was  found  that  Model  I  could  not 
adequately  account  for  the  new  experimental  results.  In  this  report,  a  new  model  is  introduced 
that  accurately  describes  the  experimental  data.  The  main  difference  between  the  models  is  the 
way  in  which  thermal  conductance  is  treated.  In  Model  I,  the  thermal  conductance  was 
approximated  by  an  effective  contact  conductance  (K=kc'A  where  kc  is  a  contact  conductivity  and 
A  is  the  cross  sectional  area).  In  the  new  model,  the  thermal  conductance  is  approximated  by  an 
effective  thermal  conductivity  (K=k  A/l  where  k  is  the  thermal  conductivity  and  1  is  a 
characteristic  length). 
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Figure  1:  The  figure  on  the  left  shows  a  1/8  in.  thick  aluminum  A1  -7075  panel  with  flat  bottom 
holes  of  1,  3/4,  and  1/2  in.  diameters  and  various  depths.  The  figure  on  the  right  shows  a  quasi¬ 
isotropic  1/8  in.  thick  graphite/epoxy  composite  with  four  flat  bottom  holes  of  1/4  in.  diameter 
and  various  depths. 


EXPERIMENTAL  PROCEDURE 

The  test  sample  was  imaged  using  a  liquid  nitrogen  cooled  Amber  Engineering  4128  InSb 
focal  plane  array  (128  x  128)  camera  with  silicon  optics  operating  in  the  3  -  5  micron  spectral 
range.  The  samples  were  thermally  excited  with  a  pair  of  xenon  arc  lamps,  each  one  powered  by 
a  5  KJoule  capacitor  bank  with  a  10  msec  discharge  time. 

The  material  used  in  this  study  was  a  1/8  in.  thick  aluminum  panel  with  various  1,  3/4, 
and  1/2  in.  diameter  flat  bottom  holes  drilled  at  depths  ranging  from  25  to  100  mil  in  steps  of 
25  mil  as  shown  in  figure  1.  The  center  to  center  distance  between  flat  bottom  holes  was  set  to  at 
least  two  diameters  to  minimize  hole  proximity  effects.  The  arc  lamps  were  positioned  to 
produce  a  uniform  distribution  of  heat  over  the  sample  surface.  The  distance  from  the  arc  lamps 
to  the  sample  was  approximately  14  in.  The  distance  from  the  camera  to  the  samples  was 
approximately  22  in.  The  data  acquisition  rate  was  43  frames/sec  and  a  total  of  100  frames  were 
collected  for  each  experiment.  Figure  1  (left)  shows  a  drawing  with  the  key  parameters  used  in 
the  study  and  model.  The  parameter  “Q”  represents  the  amount  of  energy  deposited  on  the 
surface  of  the  sample  per  unit  area.  The  parameters  ‘V’  and  “d”  represent  the  thickness  of  the 
panel  and  the  distance  from  the  surface  of  the  panel  to  the  defect  (the  words  “defect”  and  “flat 
bottom  hole”  are  used  interchangeably  in  this  report),  respectively.  The  points  “T”  (just  above  a 
defect  site)  and  “B”  (far  from  any  defect)  are  the  points  over  the  surface  of  the  sample  that  where 
used  to  calculate  the  contrast.  The  quantities  80%  and  60%  represent  the  amount  of  material 
removed  as  a  result  of  the  drilling  process.  Figure  1  (right)  shows  an  actual  frame  taken  shortly 
after  the  flash  pulse.  The  white  areas  in  the  image  photo  represent  the  flat  bottom  holes  that  tend 
to  remain  hotter  than  the  background  material  since  they  are  thinner. 
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RESULTS  AND  DISCUSSION 

Data  acquisition  was  initiated  a  few  frames  before  the  flash  discharge  in  order  to  obtain  a 
baseline  amplitude  value  for  the  sample  surface.  The  camera  frame  rate  was  set  to  43  Hz  and 
acquisition  duration  was  2.5  sec,  so  that  12-bit  digital  data  were  captured  for  the  entire  thermal 
cycle  (from  room  temperature  to  final  equilibrium  temperature).  Figure  2  shows  the  entire 
thermal  history  of  two  points  on  the  surface  of  the  panel.  The  curve  labeled  “defect”  was  taken 
from  a  point  directly  above  and  in  the  center  of  the  flat  bottom  hole  with  80%  of  material 
removed  and  it  was  referred  as  “T”  in  figure  1.  This  curve  characterizes  the  thermal  evolution  of 
a  typical  damage  site.  The  curve  labeled  “reference"  refers  to  a  point  far  away  from  any  flat 
bottom  hole  and  characterizes  the  thermal  history  of  an  undamaged  site.  Unlike  the  defect  curve, 
the  thermal  history  of  the  “reference”  point  does  not  decay  monotonically,  since  there  is  a  small 
increase  in  temperature  in  the  reference  curve  at  later  times  due  to  proximity  effects.  Both  the 
“reference”  and  “defect”  curves  exhibit  extremely  elevated  temperatures  at  early  times.  This  is  an 
artifact  produced  by  reflection  of  the  initial  flash  off  of  the  sample  and  into  the  camera.  It 
normally  disappears  when  calculating  the  contrast  curve. 


Figure  2:  (Left)  Raw  data  showing  defect  and  reference  points  and  the  resulting  contrast  curve 
(defect  -  reference).  (Right)  Typical  contrast  curves  for  various  flat  bottom  holes. 

The  difference  between  the  “reference”  and  “defect”  curves  is  termed  the  “thermal 
contrast  curve,  indicated  with  a  solid  black  line  in  figure  2  (left).  Thermal  contrast  curves  start 
and  end  with  zero  temperature  since  the  initial  and  final  equilibrium  temperatures  are  uniform 
throughout  the  entire  panel.  Figure  2  (right)  shows  contrast  curves  for  three  flat  bottom  holes 
with  1,  3/4,  and  1/2  in.  diameters  and  25  mil  from  the  surface.  It  is  clear  from  this  figure  that, 
although  the  distance  from  the  panel  surface  to  the  defect  is  the  same  for  the  three  holes,  the  peak 
contrast  temperatures  and  times  and  the  overall  shape  of  the  curves  are  significantly  different.  It 
will  be  shown  later  on  in  this  report  that  these  differences  can  be  modeled  if  lateral  heat  effects 
are  taken  into  account.  Figure  3  (left)  shows  the  peak  contrast  temperatures  for  all  flat  bottom 
holes  in  the  panel.  The  contrast  of  the  1  and  3/4  in.  diameter  holes  almost  overlap  for  all  defect 
depths.  However,  the  contrast  for  the  1/2  in.  diameter  holes  diameter  is  smaller  than  for  the  larger 
holes. 
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Figure  3:  (Left)  Peak  contrast  is  plotted  as  a  function  of  the  amount  of  material  removed.  (Right) 
The  relation  between  the  “peak  contrast”  and  the  amount  of  heat  deposited  in  the  surface  of  the 
material  is  shown  for  three  different  defects. 

An  additional  set  of  experiments  was  performed,  increased  to  study  the  relationship 
between  peak  thermal  contrast  and  amount  of  energy  delivered  to  the  sample  surface.  Figure  3 
(right)  shows  the  results  of  this  experiment  for  the  1  in.  flat  bottom  holes.  The  amount  of  energy 
deposited  on  the  surface  of  the  panel  was  controlled  by  changing  the  amount  of  charge  stored  in 
the  capacitor  banks.  Four  different  settings  were  used  which  produced  energy  levels  in  the 
following  amounts:  Q=l,  Q=2,  Q=4,  and  Q=8,  where  Q  represents  the  amount  of  energy 
deposited  per  unit  area  (the  energy  is  expressed  in  arbitrary  units).  From  figure  3  (right)  it  is  clear 
that,  for  any  given  defect,  as  the  energy  deposited  on  the  surface  increases,  the  peak  temperature 
contrast  increases  in  linear  proportion. 

CALORIMETRIC  MODEL  (ZEROTH  ORDER  APPROXIMATION) 

A  simple  theoretical  model  (zeroth  order  approximation)  was  introduced  in  our  previous 
work4.  Despite  the  fact  that  it  was  based  in  equilibrium  thermodynamics  and  that  no  lateral  heat 
effects  were  assumed,  it  correctly  accounted  for  most  of  the  observed  experimental  behavior  of 
pulsed  thermography.  Figure  4  shows  a  schematic  representation  of  the  model.  Drilled  and 
undrilled  regions  are  defined  in  the  model,  and  it  is  assumed  that  no  lateral  energy  flows  between 
these  regions.  In  the  next  section,  a  more  refined  model  that  takes  lateral  heat  transfer  effects  into 
account  will  be  derived. 


Figure  4:  This  figure  shows  a  schematic  representation  of  the  zero  lateral  flow  assumption  used 
in  the  model. 

From  simple  calorimetric  arguments,  it  can  be  established  that  qi  =  mjcTi  and  q2  =  m2cT2, 
where  qi,  mi,  and  Ti  are  the  energy  deposited  on  block  1,  the  mass  of  block  1,  and  the  final 
temperature  of  block  1  (similar  definitions  hold  for  block  2).  The  initial  temperature  of  the  panel 
can  be  assumed  to  be  0  deg.  The  mass  can  be  written  in  terms  of  the  density  as  mi  =  p  Ai  d  and 
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m2  =  p  A2  to.  If  it  is  assumed  that  the  energy  deposited  on  the  surface  of  the  panel  per  unit  area  is 
constant,  i.e.,  qi/Ai  =  q2/A2  =  Q,  then  the  final  equilibrium  temperature  difference  (or  thermal 
contrast)  between  both  blocks  Ti  -  T2  =  AT  will  be 


(\ 

u 


(1) 


This  equation  correctly  accounts  for  most  of  the  observed  experimental  behavior  of  pulsed 
thermography,  i.e., 

1 .  The  contrast  temperature  (AT)  increases  linearly  with  the  amount  of  energy  deposited 
per  unit  area  (Q). 

2.  The  higher  the  specific  heat-density  of  a  material  (pet)  the  smaller  the  contrast 
temperature  becomes  (AT  t). 

3.  The  closer  the  defect  is  to  the  surface  (d  0)  the  larger  that  the  contrast  temperature 
becomes  (AT  — >  °o). 

4.  As  the  defect  depth  approaches  the  panel  thickness  (d  t0)  the  contrast  temperature 
vanishes  (AT  — >  0). 

5.  For  a  given  defect  depth  d,  the  thicker  the  panel  (to  — >  °°)  the  larger  the  contrast 
temperature  (AT— »  Q/pcd). 

CALORIMETRIC  MODEL  (FIRST  ORDER  APPROXIMATION) 

The  previous  model  can  be  modified  to  allow  for  lateral  heat  flow  effects.  For  the  most 
general  case,  in  this  extension  of  the  model,  we  assume  that  the  in-plane  and  out-of-plane 
thermal  conductivities  are  different.  A  “poor  man’s  finite  element  approximation”  (i.e.,  three 
elements  only)  will  be  used,  but  all  the  elements  will  be  thermally  interconnected.  Figure  5 
shows  a  schematic  representation  of  the  heat  flow  model. 


Figure  5:  This  figure  shows  the  building  blocks  of  the  simple  calorimetric  model. 
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In  this  model,  it  will  be  assumed  that  all  the  energy  of  the  heat  pulse  is  absorbed  in  a  layer 
of  thickness  “d”.  This  layer  corresponds  to  the  material  above  the  defect  and  is  represented  by  the 
hatched  region  in  figure  5.  As  a  result  of  this  heating  process,  the  initial  temperature  of  this  layer 
“T0”  can  be  derived  from  simple  calorimetric  argument  from  Q  =  p  c-d-T0  where  it  is  assumed 
that  the  temperature  of  the  panel  before  the  heat  pulse  was  0  deg.  The  total  thickness  of  the  panel 
is  given  by  to  =  d  +  h  where  h  represents  the  amount  of  material  under  the  defect.  We  assume  that 
the  in-plane  material  properties  are  isotropic,  but  different  from  the  out-of-plane  properties.  The 
quantities  K  and  KL  are  the  normal  and  lateral  thermal  conductance,  respectively.  These 
quantities  can  be  expressed  in  terms  of  the  in-plane  and  out-of-plane  thermal  conductivities  of 
the  panel  as 

K=kA,/(d  +  h) 

KL  =  kLAd/R  (2) 


where  “k”  is  the  normal  (out-of-plane)  thermal  conductivity  of  the  material  (in  this  case 
aluminum),  while  “kL”  is  the  lateral  (in-plane)  thermal  conductivity  (in  the  case  of  aluminum 
k  =  kt).  These  definitions  of  the  thermal  conductances  are  the  main  differences  between  this 
model  and  the  one  introduced  previously.  Ai  and  A2  are  shown  in  figure  5  and  represent  the 
surface  area  of  the  defect  and  the  surface  area  of  rest  of  the  material,  respectively.  Ai  can  be 
written  as  Ai  =  ttR2  where  R  is  the  radius  of  the  flat  bottom  hole,  while  A2  will  be  assumed  to 
tend  to  infinity  (A2  -»  °°).  Ad  is  not  shown  explicitly  in  figure  5,  but  it  represents  the  lateral  cross 
sectional  areas  and  can  be  expressed  as  Ad  =  2rcRd.  The  set  of  differential  equations  that  define 
this  problem  are 


dT,  A,  A, 

p-A2  'd-c— =  kL  -~(T,  -T2)  +  k  ~d2  -T2) 

p-A^h-c4=k-^Th(T-T^ 


(3) 


where  Ti,  T2,  and  T2’  are  the  temperatures  of  the  different  blocks  as  shown  in  figure  5.  The 
boundary  conditions  of  the  problem  are  Ti(t=0)  =  T0 ,  T2(t=0)  =  T0 ,  T2’(t=0)  =  0.  This  set  of 
coupled  differential  equations  can  be  easily  solved  in  the  limit  when  A2  — >  00.  The  contrast  curve 
AT(t)  =  Tj(t)  -  T2(t)  obtained  from  it  is 


AT(t)  = 


_ Q _ 

pc  •  d  •  (1  -  a  +  r) 


a  k 

e"d(d+h)pcl 

\ 


1-frr  k  > 
e  d(d+h)pcl 


(4) 


where  a_  kL  Ad  d  +  h  ^  r_  £  varjat>le  “t  =  time”  should  not  be  confused  with  the 

k  Aj  R  h 

parameter  “t0  =  panel  thickness”). 

The  maximum  or  peak  thermal  contrast  can  be  calculated  by  differentiating  equation  4,  i.e., 
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which  happens  at  a  time  give  by 


AT 


_Q_ 

fi  o 

a  •  h 

pc 

{•>  tj 

.  L  . 

pc  d  t0  ,  1  +  r 

- In - 

k  1-a  +  r  a 


(5) 

(6) 


In  comparing  equations  1  and  5,  it  can  be  seen  that  the  lateral  heat  flow  effects  can  be  grouped  as 
a  multiplicative  factor  to  the  main  contrast  relation  equation  1 . 


ANALYSIS  OF  RESULTS  FOR  ALUMINUM  PANEL 


To  understand  the  effects  that  the  lateral  flow  of  heat  has  on  the  peak  thermal  contrast 
ATmax ,  equation  5  needs  to  be  studied  further.  If  it  is  assumed  that  the  material  is  isotropic,  then 
the  lateral  and  normal  thermal  conductivities  will  be  the  same.  The  parameter  “a”  can  then  be 
written  as  a  =  Ad  (d  +  h)  /  (A,R)  =  2d(t0  -  d)  /  R2 ,  and  the  lateral  heat  flow  contribution  to  the  peak 
thermal  contrast  (the  square  bracket  term  in  equation  5)  can  be  simplified  to 


1  lateral 


(d)  = 


ah 


-2—1 
ah  = 


2d(t0-d)l_Rl 


R2 


2d(t„-d) 


(7) 


Figure  6  shows  a  graph  of  equation  7  for  various  hole  diameters  as  a  function  of  defect 
depth.  From  there  it  can  be  seen  that,  for  flaws  that  are  very  close  to  the  surface  (d  ~  0)  or  for 
flaws  that  are  very  deep  in  the  material  (d  *  to),  lateral  heat  flow  effects  tend  to  disappear  in  this 
model  (or  the  lateral  heat  flow  factor  fiatmi(d)  ->  1)  and,  as  a  result,  the  entire  contribution  to  the 
peak  thermal  contrast  will  come  from  equation  1  (in  which  no  lateral  heat  flow  is  assumed).  This 
result  can  be  explained  as  follows:  In  the  limit  d  0,  the  lateral  conduction  of  heat  will  tend  to 
zero  because  the  lateral  cross  sectional  area  (27tR  d)  will  become  vanishingly  small  and  therefore 
equation  1  is  recovered.  In  the  limit  when  d  ->  to,  the  lateral  conduction  of  heat  will  tend  to  zero. 
In  this  limit,  the  temperature  gradients  will  approach  zero  and  therefore  equation  1  is  recovered. 
Finally,  when  R  «>,  the  lateral  conduction  of  heat  will  again  tend  to  zero  because  of  the  length 
over  which  the  thermal  energy  has  to  travel  is  large  and  therefore  equation  1  is  recovered. 
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Figure  6:  This  figure  shows  the  lateral  heat  flow  factor  (equation  7)  as  a  function  of  the 
depth  of  the  flaw  for  various  hole  radii. 


7 


NAWCADPAX— 99-38-TR 


Equation  5  was  used  to  fit  the  experimental  values  of  the  peak  contrast  introduced  at  the 
beginning  of  the  report.  Figure  6  shows  the  result  of  the  fits  (lines)  to  the  three  sets  of  data  points. 
Only  one  parameter  (Q/pc)  was  adjusted  to  fit  all  three  data  sets. 

Equation  4  was  used  to  fit  all  the  experimental  contrast  curves.  From  figure  7,  it  can  be 
seen  that  this  simple  model  fits  the  experimental  thermal  contrast  curves  fairly  well.  The  only 
parameter  that  was  adjusted  to  fit  all  contrast  curves  was  the  thermal  conductivity  normalized  to 
the  specific  heat-density  “k/pc”.  The  best  fit  is  for  the  data  with  60%  mass  removal.  The  model 
does  not  fit  the  time  period  around  the  maximum  for  the  80%  mass  removal  case. 


TIME  (Sec)  TIME  (Sec)  TIME  (Sec) 


Figure  7:  This  figure  shows  the  fit  (solid  lines)  to  the  experimental  contrast  curve  data  (dotted 
lines).  Each  graph  shows  the  contrast  curves  for  three  flat  bottom  holes  at  constant  depth  but  with 
different  diameter  holes. 

OUASI-ISOTROPIC  MATERIALS 

To  validate  the  model  further,  a  1/8  in.  thick  graphite  epoxy  composite  panel  was 
fabricated  with  four  0.5  in.  diameter  flat  bottom  holes  at  various  depths.  The  standard  heat  pulse 
experiment  was  performed  on  the  sample.  After  careful  analysis  of  the  data,  it  was  found  that  our 
simple  model  (equations  4,  5,  and  6)  was  not  able  to  adequately  fit  the  data.  Two  factors  are 
believed  to  be  responsible  for  the  discrepancies;  a)  the  significantly  larger  thermal  relaxation 
time  of  the  Gr/Ep  panel  compared  to  the  aluminum  panel,  and  b)  the  larger  in-plane  thermal 
conductivity  compared  to  the  out-of-plane  one.  A  new,  more  general  model  was  developed.  In 
this  model  (figure  8),  it  is  assumed  that  the  thickness  of  material  over  which  the  heat  pulse  is 
absorbed  is  a  new  variable  defined  by  “p”.  A  new  set  of  differential  equations  can  be  written  for 
this  problem  in  the  same  fashion  that  the  equations  for  the  previous  model  (equation  3)  were 
derived.  The  thermal  conductances  are  defined  in  the  same  fashion  as  before  (equation  2).  A 
closed  form  solution  can  be  found  for  the  temporal  dependence  of  the  thermal  contrast  (not 
shown  here  because  of  its  length).  A  closed  form  solution  for  the  peak  thermal  contrast  and  for 
the  time  at  which  the  peak  thermal  contrast  happens  has  not  been  found.  Figure  9  shows  the 
results  of  this  more  general  model.  Some  of  the  discrepancies  in  fit  can  be  attributed  to  the  choice 
of  position  of  the  reference  point  used  to  generate  the  background  curve. 
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Figure  8:  This  figure  shows  the  building  blocks  of  our  most  general  model. 
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Figure  9:  This  figure  shows  the  results  of  the  fit  of  the  new  model  to  the  four  flat  bottom  holes  in 
a  graphite  epoxy  composite. 
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CONCLUSIONS 

Three  simple  models  have  been  developed  that  describe  the  main  features  of  thermal 
pulse  analysis  when  applied  to  planar  flaws.  The  last  and  most  general  of  these  models  takes  into 
account  lateral  heat  conduction,  anisotropic  thermal  conductivity,  thickness,  flaw  size,  density, 
and  pulse  duration.  Equations  4, 5,  and  6  of  the  second  model  provide  a  simple,  yet  accurate, 
description  of  the  time  dependence  of  the  thermal  contrast. 
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